G Protein Activation Occurs via a Largely Universal Mechanism.

Understanding how signaling proteins like G proteins are allosterically activated is a long-standing challenge with significant biological and medical implications. Because it is difficult to directly observe such dynamic processes, much of our understanding is based on inferences from a limited number of static snapshots of relevant protein structures, mutagenesis data, and patterns of sequence conservation. Here, we use computer simulations to directly interrogate allosteric coupling in six G protein α-subunit isoforms covering all four G protein families. To analyze this data, we introduce automated methods for inferring allosteric networks from simulation data and assessing how allostery is conserved or diverged among related protein isoforms. We find that the allosteric networks in these six G protein α subunits are largely conserved and consist of two pathways, which we call pathway-I and pathway-II. This analysis predicts that pathway-I is generally dominant over pathway-II, which we experimentally corroborate by showing that mutations to pathway-I perturb nucleotide exchange more than mutations to pathway-II. In the future, insights into unique elements of each G protein family could inform the design of isoform-specific drugs. More broadly, our tools should also be useful for studying allostery in other proteins and assessing the extent to which this allostery is conserved in related proteins.

A cryptic pocket in Ebola VP35 allosterically controls RNA binding.

Protein-protein and protein-nucleic acid interactions are often considered difficult drug targets because the surfaces involved lack obvious druggable pockets. Cryptic pockets could present opportunities for targeting these interactions, but identifying and exploiting these pockets remains challenging. Here, we apply a general pipeline for identifying cryptic pockets to the interferon inhibitory domain (IID) of Ebola virus viral protein 35 (VP35). VP35 plays multiple essential roles in Ebola's replication cycle but lacks pockets that present obvious utility for drug design. Using adaptive sampling simulations and machine learning algorithms, we predict VP35 harbors a cryptic pocket that is allosterically coupled to a key dsRNA-binding interface. Thiol labeling experiments corroborate the predicted pocket and mutating the predicted allosteric network supports our model of allostery. Finally, covalent modifications that mimic drug binding allosterically disrupt dsRNA binding that is essential for immune evasion. Based on these results, we expect this pipeline will be applicable to other proteins.

SARS-CoV-2 simulations go exascale to predict dramatic spike opening and cryptic pockets across the proteome.

SARS-CoV-2 has intricate mechanisms for initiating infection, immune evasion/suppression and replication that depend on the structure and dynamics of its constituent proteins. Many protein structures have been solved, but far less is known about their relevant conformational changes. To address this challenge, over a million citizen scientists banded together through the Folding@home distributed computing project to create the first exascale computer and simulate 0.1 seconds of the viral proteome. Our adaptive sampling simulations predict dramatic opening of the apo spike complex, far beyond that seen experimentally, explaining and predicting the existence of 'cryptic' epitopes. Different spike variants modulate the probabilities of open versus closed structures, balancing receptor binding and immune evasion. We also discover dramatic conformational changes across the proteome, which reveal over 50 'cryptic' pockets that expand targeting options for the design of antivirals. All data and models are freely available online, providing a quantitative structural atlas.

SARS-CoV-2 Nsp16 activation mechanism and a cryptic pocket with pan-coronavirus antiviral potential.

Coronaviruses have caused multiple epidemics in the past two decades, in addition to the current COVID-19 pandemic that is severely damaging global health and the economy. Coronaviruses employ between 20 and 30 proteins to carry out their viral replication cycle, including infection, immune evasion, and replication. Among these, nonstructural protein 16 (Nsp16), a 2'-O-methyltransferase, plays an essential role in immune evasion. Nsp16 achieves this by mimicking its human homolog, CMTr1, which methylates mRNA to enhance translation efficiency and distinguish self from other. Unlike human CMTr1, Nsp16 requires a binding partner, Nsp10, to activate its enzymatic activity. The requirement of this binding partner presents two questions that we investigate in this manuscript. First, how does Nsp10 activate Nsp16? Although experimentally derived structures of the active Nsp16/Nsp10 complex exist, structures of inactive, monomeric Nsp16 have yet to be solved. Therefore, it is unclear how Nsp10 activates Nsp16. Using over 1 ms of molecular dynamics simulations of both Nsp16 and its complex with Nsp10, we investigate how the presence of Nsp10 shifts Nsp16's conformational ensemble to activate it. Second, guided by this activation mechanism and Markov state models, we investigate whether Nsp16 adopts inactive structures with cryptic pockets that, if targeted with a small molecule, could inhibit Nsp16 by stabilizing its inactive state. After identifying such a pocket in SARS-CoV2 Nsp16, we show that this cryptic pocket also opens in SARS-CoV1 and MERS but not in human CMTr1. Therefore, it may be possible to develop pan-coronavirus antivirals that target this cryptic pocket.

The SARS-CoV-2 nucleocapsid protein is dynamic, disordered, and phase separates with RNA.

The SARS-CoV-2 nucleocapsid (N) protein is an abundant RNA-binding protein critical for viral genome packaging, yet the molecular details that underlie this process are poorly understood. Here we combine single-molecule spectroscopy with all-atom simulations to uncover the molecular details that contribute to N protein function. N protein contains three dynamic disordered regions that house putative transiently-helical binding motifs. The two folded domains interact minimally such that full-length N protein is a flexible and multivalent RNA-binding protein. N protein also undergoes liquid-liquid phase separation when mixed with RNA, and polymer theory predicts that the same multivalent interactions that drive phase separation also engender RNA compaction. We offer a simple symmetry-breaking model that provides a plausible route through which single-genome condensation preferentially occurs over phase separation, suggesting that phase separation offers a convenient macroscopic readout of a key nanoscopic interaction.

SARS-CoV-2 Nsp16 activation mechanism and a cryptic pocket with pan-coronavirus antiviral potential.

Coronaviruses have caused multiple epidemics in the past two decades, in addition to the current COVID-19 pandemic that is severely damaging global health and the economy. Coronaviruses employ between twenty and thirty proteins to carry out their viral replication cycle including infection, immune evasion, and replication. Among these, nonstructural protein 16 (Nsp16), a 2'-O-methyltransferase, plays an essential role in immune evasion. Nsp16 achieves this by mimicking its human homolog, CMTr1, which methylates mRNA to enhance translation efficiency and distinguish self from other. Unlike human CMTr1, Nsp16 requires a binding partner, Nsp10, to activate its enzymatic activity. The requirement of this binding partner presents two questions that we investigate in this manuscript. First, how does Nsp10 activate Nsp16? While experimentally-derived structures of the active Nsp16/Nsp10 complex exist, structures of inactive, monomeric Nsp16 have yet to be solved. Therefore, it is unclear how Nsp10 activates Nsp16. Using over one millisecond of molecular dynamics simulations of both Nsp16 and its complex with Nsp10, we investigate how the presence of Nsp10 shifts Nsp16's conformational ensemble in order to activate it. Second, guided by this activation mechanism and Markov state models (MSMs), we investigate if Nsp16 adopts inactive structures with cryptic pockets that, if targeted with a small molecule, could inhibit Nsp16 by stabilizing its inactive state. After identifying such a pocket in SARS-CoV-2 Nsp16, we show that this cryptic pocket also opens in SARS-CoV-1 and MERS, but not in human CMTr1. Therefore, it may be possible to develop pan-coronavirus antivirals that target this cryptic pocket. STATEMENT OF SIGNIFICANCE: Coronaviruses are a major threat to human health. These viruses employ molecular machines, called proteins, to infect host cells and replicate. Characterizing the structure and dynamics of these proteins could provide a basis for designing small molecule antivirals. In this work, we use computer simulations to understand the moving parts of an essential SARS-CoV-2 protein, understand how a binding partner turns it on and off, and identify a novel pocket that antivirals could target to shut this protein off. The pocket is also present in other coronaviruses but not in the related human protein, so it could be a valuable target for pan-coronavirus antivirals.

Mechanism of Nucleotidyltransfer Reaction and Role of Mg2+ Ion in Sugar Nucleotidyltransferases.

Sugar nucleotidyl transferases (SNTs) catalyze nucleotidyltransfer reactions to form sugar-nucleotides and pyrophosphate in the presence of two Mg2+ ions (Mg2+A and Mg2+B). We unveil the mechanism and free energetics of nucleotidyl transfer reaction in an SNT called GlmU through hybrid quantum mechanics-molecular mechanics molecular dynamics simulations and free energy calculations. The study identifies the roles of the active site residues and the Mg2+ ions in catalyzing the reaction. Of great significance, we are able to compare the free energy barrier for the reaction with that for the Mg2+-assisted release of the product (i.e., pyrophosphate) into the solution, shedding light on the general mechanistic and kinetic aspects of catalysis by SNTs.

SARS-CoV-2 Simulations Go Exascale to Capture Spike Opening and Reveal Cryptic Pockets Across the Proteome.

SARS-CoV-2 has intricate mechanisms for initiating infection, immune evasion/suppression, and replication, which depend on the structure and dynamics of its constituent proteins. Many protein structures have been solved, but far less is known about their relevant conformational changes. To address this challenge, over a million citizen scientists banded together through the Folding@home distributed computing project to create the first exascale computer and simulate an unprecedented 0.1 seconds of the viral proteome. Our simulations capture dramatic opening of the apo Spike complex, far beyond that seen experimentally, which explains and successfully predicts the existence of 'cryptic' epitopes. Different Spike homologues modulate the probabilities of open versus closed structures, balancing receptor binding and immune evasion. We also observe dramatic conformational changes across the proteome, which reveal over 50 'cryptic' pockets that expand targeting options for the design of antivirals. All data and models are freely available online, providing a quantitative structural atlas.

A revised mechanism for (p)ppGpp synthesis by Rel proteins: The critical role of the 2'-OH of GTP.

Bacterial Rel proteins synthesize hyperphosphorylated guanosine nucleotides, denoted as (p)ppGpp, which by inhibiting energy requiring molecular pathways help bacteria to overcome the depletion of nutrients in its surroundings. (p)ppGpp synthesis by Rel involves transferring a pyrophosphate from ATP to the oxygen of 3'-OH of GTP/GDP. Initially, a conserved glutamate at the active site was believed to generate the nucleophile necessary to accomplish the reaction. Later this role was alluded to a Mg2+ ion. However, no study has unequivocally established a catalytic mechanism for (p)ppGpp synthesis. Here we present a revised mechanism, wherein for the first time we explore a role for 2'-OH of GTP and show how it is important in generating the nucleophile. Through a careful comparison of substrate-bound structures of Rel, we illustrate that the active site does not discriminate GTP from dGTP, for a substrate. Using biochemical studies, we demonstrate that both GTP and dGTP bind to Rel, but only GTP (but not dGTP) can form the product. Reactions performed using GTP analogs substituted with different chemical moieties at the 2' position suggest a clear role for 2'-OH in catalysis by providing an indispensable hydrogen bond; preliminary computational analysis further supports this view. This study elucidating a catalytic role for 2'-OH of GTP in (p)ppGpp synthesis allows us to propose different mechanistic possibilities by which it generates the nucleophile for the synthesis reaction. This study underscores the selection of ribose nucleotides as second messengers and finds its roots in the old RNA world hypothesis.

The SARS-CoV-2 nucleocapsid protein is dynamic, disordered, and phase separates with RNA.

The SARS-CoV-2 nucleocapsid (N) protein is an abundant RNA binding protein critical for viral genome packaging, yet the molecular details that underlie this process are poorly understood. Here we combine single-molecule spectroscopy with all-atom simulations to uncover the molecular details that contribute to N protein function. N protein contains three dynamic disordered regions that house putative transiently-helical binding motifs. The two folded domains interact minimally such that full-length N protein is a flexible and multivalent RNA binding protein. N protein also undergoes liquid-liquid phase separation when mixed with RNA, and polymer theory predicts that the same multivalent interactions that drive phase separation also engender RNA compaction. We offer a simple symmetry-breaking model that provides a plausible route through which single-genome condensation preferentially occurs over phase separation, suggesting that phase separation offers a convenient macroscopic readout of a key nanoscopic interaction.

Mechanism of Mg2+-Accompanied Product Release in Sugar Nucleotidyltransferases.

The nucleotidyl transfer reaction, catalyzed by sugar nucleotidyltransferases (SNTs), is assisted by two active site Mg2+ ions. While studying this reaction using X-ray crystallography, we captured snapshots of the pyrophosphate (product) as it exits along a pocket. Surprisingly, one of the active site Mg2+ ions remains coordinated to the exiting pyrophosphate. This hints at the participation of Mg2+ in the process of product release, besides its role in catalyzing nucleotidyl transfer. These observations are further supported by enhanced sampling molecular dynamics simulations. Free energy computations suggest that the product release is likely to be rate limiting in SNTs, and the origin of the high free energy barrier for product release could be traced back to the "slow" conformational change of an Arg residue at the exit end of the pocket. These results establish a dual role for Mg2+, and propose a general mechanism of product release during the nucleotidyl transfer by SNTs.

A salt-bridge stabilized C-terminal hook is critical for the dimerization of a Bowman Birk inhibitor.

Legume Bowman-Birk inhibitors (BBIs) that inhibit mammalian proteases exist as dimers in solution. The structural basis governing dimerization of HGI-III (horsegram seed BBI) was investigated. An intra-monomer salt bridge (D76-K71) stabilizes an atypical hook-like conformation at the C-terminus. We postulate that this hook, positions D75 to enable an inter-monomer salt-bridge D75(a)-K24(b), which results in dimerization. We verify this by K71A and D76A mutations of HGI-III. The mutants were both monomers, likely due to destabilization of the C-terminal hook. Dimerization was sustained in a double mutant K71D/D76K that was anticipated to form a similar hook critical for dimerization. Conversely, K24(b) that interacts with D75(a) of the loop is the specificity determining residue that interacts with trypsin to inhibit its activity. The inter-monomer salt bridge D75(a)-K24(b) must be disrupted for the inhibition of trypsin, requiring HGI-III to transition into a monomer. Size exclusion studies and a model of HGI-III-trypsin complex support this notion. Interestingly, isoforms of the inhibitor present in germinated seeds (HGGIs) are monomers; and most strikingly, the C-termini of these inhibitors are truncated with the loss the C-terminal hook critical for dimerization. The tendency of HGI-III to self-associate seems to relate to its physiological function of a storage protein.

GlmU (N-acetylglucosamine-1-phosphate uridyltransferase) bound to three magnesium ions and ATP at the active site.

N-Acetylglucosamine-1-phosphate uridyltransferase (GlmU), a bifunctional enzyme exclusive to prokaryotes, belongs to the family of sugar nucleotidyltransferases (SNTs). The enzyme binds GlcNAc-1-P and UTP, and catalyzes a uridyltransfer reaction to synthesize UDP-GlcNAc, an important precursor for cell-wall biosynthesis. As many SNTs are known to utilize a broad range of substrates, substrate specificity in GlmU was probed using biochemical and structural studies. The enzymatic assays reported here demonstrate that GlmU is specific for its natural substrates UTP and GlcNAc-1-P. The crystal structure of GlmU bound to ATP and GlcNAc-1-P provides molecular details for the inability of the enzyme to utilize ATP for the nucleotidyltransfer reaction. ATP binding results in an inactive pre-catalytic enzyme-substrate complex, where it adopts an unusual conformation such that the reaction cannot be catalyzed; here, ATP is shown to be bound together with three Mg2+ ions. Overall, this structure represents the binding of an inhibitory molecule at the active site and can potentially be used to develop new inhibitors of the enzyme. Further, similar to DNA/RNA polymerases, GlmU was recently recognized to utilize two metal ions, MgA2+ and MgB2+, to catalyze the uridyltransfer reaction. Interestingly, displacement of MgB2+ from its usual catalytically competent position, as noted in the crystal structure of RNA polymerase in an inactive state, was considered to be a key factor inhibiting the reaction. Surprisingly, in the current structure of GlmU MgB2+ is similarly displaced; this raises the possibility that an analogous inhibitory mechanism may be operative in GlmU.

Identification and characterization of a hitherto unknown nucleotide-binding domain and an intricate interdomain regulation in HflX-a ribosome binding GTPase.

A role for HflX in 50S-biogenesis was suggested based on its similarity to other GTPases involved in this process. It possesses a G-domain, flanked by uncharacterized N- and C-terminal domains. Intriguingly, Escherichia coli HflX was shown to hydrolyze both GTP and adenosine triphosphate (ATP), and it was unclear whether G-domain alone would explain ATP hydrolysis too. Here, based on structural bioinformatics analysis, we suspected the possible existence of an additional nucleotide-binding domain (ND1) at the N-terminus. Biochemical studies affirm that this domain is capable of hydrolyzing ATP and GTP. Surprisingly, not only ND1 but also the G-domain (ND2) can hydrolyze GTP and ATP too. Further; we recognize that ND1 and ND2 influence each other's hydrolysis activities via two salt bridges, i.e. E29-R257 and Q28-N207. It appears that the salt bridges are important in clamping the two NTPase domains together; disrupting these unfastens ND1 and ND2 and invokes domain movements. Kinetic studies suggest an important but complex regulation of the hydrolysis activities of ND1 and ND2. Overall, we identify, two separate nucleotide-binding domains possessing both ATP and GTP hydrolysis activities, coupled with an intricate inter-domain regulation for Escherichia coli HflX.

Crystal structures identify an atypical two-metal-ion mechanism for uridyltransfer in GlmU: its significance to sugar nucleotidyl transferases.

N-Acetylglucosamine-1-phosphate uridyltransferase (GlmU), exclusive to prokaryotes, is a bifunctional enzyme that synthesizes UDP-GlcNAc-an important component of the cell wall of many microorganisms. Uridyltransfer, one of the reactions it catalyzes, involves binding GlcNAc-1-P, UTP and Mg(2+) ions; however, whether one or two ions catalyze this reaction remains ambiguous. Here, we resolve this using biochemical and crystallographic studies on GlmU from Mycobacterium tuberculosis (GlmU(Mtb)) and identify a two-metal-ion mechanism (mechanism-B). In contrast to well-established two-metal mechanism (mechanism-A) for enzymes acting on nucleic acids, mechanism-B is distinct in the way the two Mg(2+) ions (Mg(2+)A and Mg(2+)B) are positioned and stabilized. Further, attempts to delineate the roles of the metal ions in substrate stabilization, nucleophile activation and transition-state stabilization are presented. Interestingly, a detailed analysis of the available structures of sugar nucleotidyl transferases (SNTs) suggests that they too would utilize mechanism-B rather than mechanism-A. Based on this, SNTs could be classified into Group-I, which employs the two-metal mechanism-B as in GlmU, and Group-II that employs a variant one-metal mechanism-B, wherein the role of Mg(2+)A is substituted by a conserved lysine. Strikingly, eukaryotic SNTs appear confined to Group-II. Recognizing these differences may be important in the design of selective inhibitors against microbial nucleotidyl transferases.

Substrate-bound crystal structures reveal features unique to Mycobacterium tuberculosis N-acetyl-glucosamine 1-phosphate uridyltransferase and a catalytic mechanism for acetyl transfer.

N-acetyl-glucosamine-1-phosphate uridyltransferase (GlmU), a bifunctional enzyme involved in bacterial cell wall synthesis is exclusive to prokaryotes. GlmU, now recognized as a promising target to develop new antibacterial drugs, catalyzes two key reactions: acetyl transfer and uridyl transfer at two independent domains. Hitherto, we identified GlmU from Mycobacterium tuberculosis (GlmU(Mtb)) to be unique in possessing a 30-residue extension at the C terminus. Here, we present the crystal structures of GlmU(Mtb) in complex with substrates/products bound at the acetyltransferase active site. Analysis of these and mutational data, allow us to infer a catalytic mechanism operative in GlmU(Mtb). In this S(N)2 reaction, His-374 and Asn-397 act as catalytic residues by enhancing the nucleophilicity of the attacking amino group of glucosamine 1-phosphate. Ser-416 and Trp-460 provide important interactions for substrate binding. A short helix at the C-terminal extension uniquely found in mycobacterial GlmU provides the highly conserved Trp-460 for substrate binding. Importantly, the structures reveal an uncommon mode of acetyl-CoA binding in GlmU(Mtb); we term this the U conformation, which is distinct from the L conformation seen in the available non-mycobacterial GlmU structures. Residues, likely determining U/L conformation, were identified, and their importance was evaluated. In addition, we identified that the primary site for PknB-mediated phosphorylation is Thr-418, near the acetyltransferase active site. Down-regulation of acetyltransferase activity upon Thr-418 phosphorylation is rationalized by the structures presented here. Overall, this work provides an insight into substrate recognition, catalytic mechanism for acetyl transfer, and features unique to GlmU(Mtb), which may be exploited for the development of inhibitors specific to GlmU.